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Abstract 

We study the response of an infinite, asymptotically static M = 4 plasma to a generic 
localized source in the probe approximation. At large distances, the energy momentum 
tensor of the plasma includes a term which satisfies the constitutive relations of linearized 
hydrodynamics, but it can also include a non-hydrodynamical term which contributes at the 
same order as viscous corrections, or even at leading order in some cases. The conditions for 
the appearance of a laminar wake far behind the source and its relevance for phenomeno- 
logical models used to explain di-hadron correlations are discussed. We also consider the 
energy momentum tensor near the source, where the hydrodynamical approximation can 
be expected to break down. Our analysis encompasses a wide range of sources which are 
localized in the bulk of AdS, including trailing strings, mesonic and baryonic configurations 
of strings, and point particles. 
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1 Motivation 

As a parton moves through a quark-gluon plasma, it loses energy to the medium. It is inter- 
esting to ask where this energy goes. At sufficiently large length scales, the main channels 
available are hydrodynamical. There are two different hydrodynamical modes. One is sound. 
Assuming the parton is moving at a speed v greater than the speed of sound c s , the signature 
of energy loss into the sound mode is a sonic boom. The other mode is dispersive, having 
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to do with the formation of laminar wakes (or diffusion wakes — we will not attempt to draw 
a distinction between them). Neglecting effects related to the expansion or other collective 
flows of the medium, the laminar wake is a stream of fluid behind the parton, moving in 
the same direction as the parton. Once the plasma hadronizes, we can expect that a sonic 
boom will lead to enhanced particle production at the Mach angle 9 = cos~ l c s /v, while a 
hadronized wake will lead to intensified particle production in the direction of the parton's 
motion [lj [2] . 

Evidence for such effects may be obtained from histograms of the azimuthal angle between 
pairs of energetic hadrons produced in heavy ion collisions OH]. In [3] it was observed that 
for certain momenta the two-point correlation function between jets emitted from the plasma 
is peaked at an azimuthal angle of roughly 7r± 1.2 radians and has a minimum at ir radians. 
This is described as "jet-splitting" . Jet-splitting is suggestive of a sonic boom: a standard 
interpretation is that one energetic hadron (the "trigger" or "near-side" hadron) came from 
a hard parton that exited the plasma without losing much energy, and the other one (the 
"associated" or "away-side" hadron) was produced from the sonic boom caused by another 
hard parton whose momentum was opposite the first, at least in azimuthal angle. In [I], 
with more inclusive momentum cuts (and also the greater rapidity acceptance characteristic 
of STAR), instead of jet-splitting, a broad peak was found, centered around ir\^ 

It is challenging to make an unambiguous connection between these results and hydro- 
dynamics, but there are several notable efforts. For example, in [2] it was shown, using a 
hydrodynamical model and Cooper-Frye hadronization, that jet-splitting doesn't occur un- 
less the diffusion wake is suppressed relative to the sonic boom; and in [BJ, jet-splitting, in 
approximate agreement with PHENIX di-hadron correlators, was predicted using a model 
where three quarters of the energy goes into the sonic boom. The discussion so far is any- 
thing but an exhaustive account of either the experimental or the theoretical literature on 
medium-induced modifications of jet structure. Recent brief discussions can be found in 
[3, [8] , while a broader treatment with more extensive references is included in [9] . 

In both [2] and in [6J the authors have tuned the relative amount of energy going into 
sound modes and diffusion modes by hand. Indeed, it is challenging to predict from QCD 
the relative strength of the sonic boom and the diffusion wake produced by a hard parton — 
assuming, of course, that sonic booms and diffusion wakes are the right language for describ- 

*A recent analysis [5] of STAR data using three-point functions exhibits a triple peak structure, where 
a central peak at A0i = A<p2 = ir is accompanied by two peaks of equal heights at A0x ~ 7r ± 1.4 and 
A(f>2 ~ 7r =p 1.4. These latter two peaks are consistent with a sonic boom, while the central peak may be 
evidence for a diffusion wake. 



2 



ing the energy loss at length scales significantly larger than 1 fm. The essential difficulty is 
in connecting the short- distance physics, which is perturbative (or at least partly perturba- 
tive) and the long-distance hydrodynamical regime. It might therefore be enlightening to 
consider similar phenomena in M = 4 super- Yang- Mills theory, where the methods of the 
gauge-string duality [10l [11], [12] allow detailed calculations which become reliable in the large 
N, large Oym-N limit. Indeed, a number of papers [131 dH d51 ESI CEO EES! EES] are devoted to 
studying the sonic boom and diffusion wake produced by the trailing string of [2Ql 121] , which 
represents a heavy quark moving at constant velocity through a thermal medium of M = 4 
gauge theory. A strong diffusion wake is predicted in these works, similar to "scenario 1" 
of [2], which led to no jet-splitting after Cooper- Frye hadronization. On the other hand, in 
[22] . we showed that a string configuration representing a heavy-quark meson as described 
in [231 121] produces a sonic boom but no diffusion wake. 

The aim of this paper is to get a stronger foothold on the hydrodynamic behavior of the 
M = 4 plasma due to a generic probe source, with an eye toward the phenomenologically 
interesting question of when a diffusion wake arises and how strong it is relative to the sonic 
boom. Our main result, which characterizes the linear response of the stress-tensor of the 
boundary theory to the source and provides a comparison to linear hydrodynamics, is pre- 
sented in section [2J Hydrodynamic modes in the strongly coupled J\f = 4 plasma have been 
widely studied following [25] , and recent works include [261 [23 CEEl HSl [281 [22] , where more 
extensive references may be found. In section [3] and H] we give details of our computation. In 
section [5] we discuss the near field of the source, which seems non-hydrodynamical in nature 
[5U] . In section [H] we apply our methods to strings and point particle sources. 



2 Summary of results 



We start with an action 



S = -L f d 5 x y/^G 
2 «i J 
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where Sm describes a localized source. Without the source, one finds the standard transla- 
tionally invariant AdiS'5-Schwarzschild solution to the Einstein equations 



dz 2 



ds 2 = a 2 -hdt 2 + ^ dx "i + —T 



(2) 



i=l 
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where 



L 



z 



1 



«(«) = -, /i(z) = l- - . (3) 
2; \ z , 



The asymptotically AdS boundary is located at z = and the black hole horizon is at z = z . 
The AdS/CFT duality relates strings moving in an asymptotically AdSsxS 5 background to 
M = 4 Super Yang Mills theory. We have omitted the S 5 part of the metric, assuming 
that the equations of motion of the source and the Einstein equations can be consistently 
truncated to the AdS5 directions. In [31] the role of the motion of a string-source along the 
S 5 directions was studied in a rotating black hole background. 

Treating a source in the probe approximation means determining its motion by extrem- 
izing Sm in the background (j2J) and then solving the linearized Einstein equations in the 
presence of the source to find a perturbed metric. Let the energy momentum tensor of the 
source be denoted J^(t,x, z), where indices /1 and v run over the five dimensions of AdS 5- 
Schwarzschild. The tensor J^ u is defined by the first variation of Sm with respect to the 
metric: 



^3 J d'xV^GdG^J^, (4) 
and the Einstein equations (before linearization) read 



SS M 



p Ipr* —r —^JLi (5) 

The factors of 1/L 3 in front of J Mi , in and (jSJ) are included for later convenience. 

The gauge-string duality relates the metric perturbations near the boundary of AdS^ to 
the expectation value of the energy-momentum tensor of the boundary gauge theory, which 
we denote by (T mn ). The indices m and n run from to 3, and 5 represents the radial AdS 
direction: that is, x 5 = z. Our first main result is a generalization of the observation, first 
made quantitatively in [13J, that T mn is not necessarily conserved in the presence of a source: 
instead, 

d m (T mn (t,x)) = J^(t,x), (6) 
where Jjy(t,x) is defined through 

00 

J^(t, x,z) = Yl J ^M)2 a > (7) 

a=ao 

where a is integer- valued and ao > — 1. There are some caveats to (J6j): if Jg 5 7^ then 
(T mn ) will have a trace anomaly. In addition, (T mn ) can have divergent terms which can 
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contribute to the non- conservation law ([6]). These divergences can usually be ignored in 
phenomenological applications because their contribution to (T mn ) is confined to the location 
of the probe. In section [3] and appendix [C] we derive the result §6§ and explain this caveat 
in more detail. 

Perhaps the most phenomenologically interesting information about the response of the 
medium to the source comes from the long distance asymptotics, where (T mn ) is expected to 
approximately satisfy the constitutive .Miens of hydrodyna^a Our second main result 
characterizes the large distance behavior of (T mn ) in Fourier space, 

(Tmn(t, 2)> = J ^Wf e ~ Mt3{Tmn{U ' ^ ' (8) 

Still within the probe approximation, we find that the small k behavior of the boundary 
theory stress-energy tensor (T mn ) is given by 

(-^mn) {Tmn) bath 

+ T% n dl ° + T mn + (corrections) (9) 

where 

7T 2 

(T mn ,) bath = — (N 2 - 1)T 4 diag (3, 1,1,1) (10) 

o 

is the energy momentum of the thermal bath, N is the number of colors in the Af = 4 gauge 
theory, and T is the temperature. The correction terms, labeled (corrections) in (Q, are 
smaller than T m ydro by a factor of at least 0(k). If a multipole expansion of the drag force 
density, J^, includes a monopole term then these corrections are smaller than T m y r f vo by a 
factor of 0(k 2 ). The tensor jF mn appearing in the third term of Q is defined as follows: 



T 



c 



3%)(0 



•7~m0 — ^Om — , (11' 



where F [f(z)] means the finite part of f(z) when taking the z — > limit, and the index i 
runs from 1 to 3. The notation (ij) means to take the symmetric traceless part of a 3 x 3 
matrix: 

X (ij) = ^i Xij + Xji ^ ~ 3 S V S9kX 9h ■ ( 12 ) 

As in ([6]), we have discarded from (Q divergent contributions to the stress-energy tensor, a 



2 An interesting recent study [30j explores the degree to which the non-linear constitutive relations are 
maintained in the near-field regime of the trailing string, building upon results of [321 1141 [33] . Comparisons 
were also made for the trailing string between string theory results and linearized hydrodynamics in [16j . 
and between string theory and a regulated version of linearized hydro results in |15j . 
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trace anomaly term and additional finite terms which are related to logarithmic divergences. 
We discuss these terms in detail in sections EJ @] and appendix O where we holographically 
renormalize the stress-energy tensor. 

The tensor T^ dro appearing in the second term of (Q is defined by two requirements. 
First, it should satisfy the constitutive relations of linearized hydrodynamics, which means 
that it can be written in the form 

S l \ 

1 . 1 Q O ' ( 13 ) 

/ 

where we have again used the notation (|T2|) . Note that e is the deviation of (T 00 ) from 
its value in the thermal bath, whereas Si is the complete Poynting vector, simply because 
neither the bath nor T mn contributes to it. The second requirement on T^ TO is that it 
should be traceless and should obey the conservation equation that follows from plugging 
P) into ©: that is, 

gm T Wro = F hydro = _ gm ^ (14) 

The explicit values of e and S after implementing conditions (1131) and (1141) can be found in 
( 1631) . We could summarize by saying that T^ vo captures most of the leading asymptotics 
at small k, and T mn characterizes the deviations from the leading hydro dynamical behavior. 
Deviations of the quark gluon plasma from hydrodynamical behavior was also observed in 
[2H] when the gauge theory was placed on a sphere. We note that the deviations we observe 
here are different in nature from those in since they are associated with excitations 
created by the probe source and not with inherent properties of the fluid. In [TTJ, [22j [16] 
these deviations have been studied for the specific case of a heavy quark and heavy meson 
source. We emphasize that (1131) has been obtained by solving the equations of motion for 
the metric perturbations in a black hole background — we have not assumed its form a priori. 
Thus, the constitutive relations we've obtained in (fl~3|) provide an alternate derivation of the 
standard result t]/s = I/Att [25]. 

Obtaining our two main results, and requires some technical assumptions. To 
obtain (jH]), we must assume that ([7]), a power series expansion of near the boundary, 
exists. To obtain (Q, we must additionally assume that a near boundary power series 



rpmn 
hydro 
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expansion exists also for the Fourier coefficients of J^ u {oj, k, z) defined through 

J,u(t, x ,z) = J ^ ^5 k, z) , (f 5) 

and that the source is in some sense localized in the x m directions. The latter assumption 
may be made precise in one of three ways: 

1. Either Jfj, v (uj, k, z) should be analytic in to and k\ 

2. Or it should be expressible in the form 

J^iuj, k, z) = A(u, k)j^(u, k, z) , (16) 

where j^iyj, k, z) is analytic in uj and k but A (a;, k) need not be; 

3. Or it should be expressible as a sum of terms of the form shown on the right hand side 
of (HSJ). 

In all three cases one must require that J M ^(t, x, z) is conserved: that is, D^J^ = with 
the covariant derivative in the background ([2]). We must also make some mild assumptions 
about the behavior of near the horizon to ensure that a non-singular metric response 
there is possible. The precise form of these assumptions will be discussed in section I3TT1 and 
H 

3 The boundary stress tensor and its non-conservation 

Our overall goal is to evaluate the energy momentum tensor of the boundary theory (T mn ) in 
terms of the energy momentum tensor of a probe-source moving in the ArfSs-Schwarzschild 
background. The focus of this section is to extract the maximum possible information from 
a near-boundary analysis and to put it in a form that will make it easy to consider the large 
distance asymptotics, which we do in section HI In section [3TT1 we discuss the consequences of 
the conservation equation for the energy momentum tensor J^ v of the probe. In section 13.21 
we perform a near-boundary analysis of the perturbations of the five-dimensional metric and 
explain how (T mn ) can be extracted from them. In section [3731 we use constraints on (T mn ) 
to give a particular parametrization of it in terms of five undetermined coefficients, whose 
subsequent study will occupy us in section |U In section \5\ we study the near field of the 
source using a different parametrization from the one given in 13.31 
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3.1 Constraints on the source terms 

Consider a probe-source moving in the AdSs-Schwarzschild background (j2J). We require that 
the stress tensor of the probe-source J^ u is conserved: 

iV^O, (17) 

where is the usual covariant derivative in curved spacetime. If ffTTl) fails, the Einstein 
equations will be inconsistent. One way to see this inconsistency is to pass to an axial gauge 
where five of the fifteen independent components of the metric are set to zero. There are 
nevertheless fifteen Einstein equations for the remaining ten components of the metric, and 
five of them are first order constraints whose consistency with the ten second order equa- 
tions of motion depends on (IT7|) being true. In terms of the Fourier coefficients J Ml /(o;, k, z) 
introduced in (TT5|) . the conservation equation (1171) takes the form 

iUJ T i*7 t {Jm5a 3 h)' 

T^mO ~r IKiJmi — o {loa.) 

h cr 

™ 7 , it j T a> t w ( j ^ a2h3/2 y n«w 



where we have defined 



(-W h k 2 fc 3 ) (19) 



and used the labels m = 0, 1, 2, 3 and i = l,2, 3. 

Consider a near boundary series expansion of the Fourier coefficients Jfj, v (u),k,z) in 
similar to (JTj), 

oo 

J, u (u,k,z) = Y J 4»^^)z a . (20) 

a=ao 

We will assume that a runs over integers. Near the boundary ~ l/^ 2 , so we must 
take a > —2 in order for the contribution of J^(uj, k, z) to the Einstein equations to be 
less divergent near z = than the cosmological constant term j^G^. If < —2, then 
generically the spacetime will no longer be asymptotically anti-de Sitter. Plugging (J7j) into 
the conservation equations ffTS]) . one finds that = 0, and that the first few coefficients 



3 In section [2] we claimed that the existence of a series expansion of the Fourier coefficients in small z is 
not required to obtain ([6]). This is correct as can be seen by exchanging ik m with d m in the analysis leading 
to (|34|) and (|35|) . Nevertheless, we have chosen to work in Fourier space so that the results we obtain here 
may be easily referred to from later sections where we do require a series expansion in small z for the Fourier 
coefficients of J^ v . 



S 



JjS are related by the equations 

Jnm ' = (3 — a )Jn5 

4r 1} (3 - a) = + ifc m J,Sr 2) (21) 

where a = 0, ... ,3. The components J55 and remain undetermined by these relations. 
They are integration constants of (ITTj) . In (I2ip and in the rest of this work, the index m is 
raised using the Minkowski metric. 



3.2 Perturbations of the five-dimensional metric 

In order to compute (T mn ) using the gauge-string duality we need to find the small z asymp- 
totics of the linear response of the metric to the source. We write the perturbations as 
follows: 

2 2 

n _ ^yAdSBH , K 5 a tt (oi\ 
^fMU — tr^v + ~2lfi ^ ' ^ ' 

where G^ SBH dx fJ 'dx u is the AdSs-Schwarzschild line element (j2J), and is small compared 
to G^ SBH . The normalization factors in the second term in ( 1221) are chosen so that if 
= H$z 4 + 0(z 5 ) and we work in a gauge where = 0, then the holographic stress 
tensor is 

(Tmn) = (T mn )ba.th + ■ (^3) 

To determine H^ u , one must solve the linearized version of the Einstein equations (|SJ), 
which we will write formally as 

V^ pa H^ = . (24) 

The equations where \i and v run from to 3 comprise ten second order equations of motion, 
and the equations with /i = 5 give five first order constraint equations. The most general 
solution to the second order equations of motion involves twenty integration constants. Ten 
set the value of H mn (0) and are related to deformations of the boundary theory metric. We 
insist that this metric should not be deformed, and so H mn (0) = 0. If we use the following 
expansion (c.f. [34"]): 

H mn (z) = H^ n z + H^ n z + H^ n z + H^ n z In z/L + H^z + . . . , (25) 

then the other ten integration constants are given by Hmn- The Hmn with a < 4, and also 
Hmn, may be determined by solving the second order equations of motion perturbatively in 
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z. We find for a = 1, 2 or 3, 

= - fj( a-2) --77 J {a ~ 2)s - — ( R {a ~ 2) - -R {a - 2) n X) (26a) 

11 mn a(4 — a) \ n 3 1 k 2 I nm g 'mm l y ; ^ Ud i 

and also 

#(4) _ _ ( j(2) _ 1 j(2) S _^fn(2)_lr,(2) ^ ( 2 6b) 

where i? nm and i? are the Ricci Scalar and Ricci tensor which follow from the metric g mn = 

2 

Vmn + jiflH mn expanded to linear order in H mn . Following the notation in (0) and ([2"5|) . we 
have defined 

J R = i? (1) 2 + J R (2) 2 2 + .... (28) 

So, for example, Rmn is a linear combination of second derivatives of Hmn. 

Once the H mn are known, we can compute the boundary energy momentum tensor by 
varying the action with respect to the boundary value of the metric. We carry out this 
procedure in appendix O The resulting stress tensor takes the form 

{Tmn) — (Tmn/bstix + V^mn) + (^nrn) 5 (29) 

where (5T mn ) is the finite part of the fluctuation of the stress-tensor above the background 
value of the plasma and (T mn ) is a divergent contribution to the stress energy tensor, whose 
explicit form is given in (11601) . Usually, divergent terms in the stress energy tensor can be 
removed by introducing appropriate counter-terms via holographic renormalization. The 
infinities which appear in fl29l) are of a different nature: we expect that they are associated 
with the probe source having parameters which formally diverge. For example, string con- 
figurations which end on the asymptotically AdS boundary are dual to infinitely massive 
quarks. We show in section [6] that all the divergent terms in the boundary theory stress 
tensor (T mn ) are a result of the infinite mass of the quark. More generally, to the extent that 
the J mn are localized in space, the divergent terms in (T mri ) are likewise localized, meaning 
that we may treat them as divergent contact terms supported at the location of the probe. 
See appendix [C] for an extended discussion of these divergences. 

The non-divergent part of the fluctuation of the stress-energy tensor above the back- 
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ground value of the plasma is given by 

3.3 Constraints on the boundary stress tensor 

To evaluate (ST mn ), we need to calculate the Hmn's in ( 1251) . Five of the ten coefficients Hmn 
can be determined by perturbatively solving the five first order constraint equations. From 
the 55 equation, 

rrm (4) _ 1 jm (2) _ 1 t(2) /qin 

and from the m5 equations, 

zk m H^ = jS +i k n H^, (32) 

where we made use of ( 1211) and (]26p . Had the conservation condition fTiTI) not been satisfied, 
the constraint equations would have been inconsistent with the second order equations of 
motion. In light of the constraints (13T!) and (|32|) . there are only five independent components 
of Hmn- They are determined by solving the Einstein equations with inf ailing boundary 
conditions at the horizon. We will explain how to do this explicitly in section H] in the 
context of a small k expansion (and in section [5] for a large k expansion.) Here we wish to 
avoid any use of an expansion in small or large k, but we will solve the equations fl3T|) and fl32l) 
in a way that will make it as easy as possible to study a small k expansion later. A good 
strategy is to choose the undetermined components to be the traceless spatial part of H mn , 
which we denote as H^y The reason this is a good strategy is that when hydrodynamics 
is valid, the components (5T(jj)) of the stress tensor are usually one order in k smaller than 
the other components for small k: they arise from viscous effects in the regime of linearized 
hydrodynamics. It is straightforward to use (13"T|) . (13"2"1) . and 

^!d(2) _ _ r(2)m jk m —RW - -hk T^ m (33^1 

together with fl3"Ul) . to show that 

ik m (5T mn ) = jf 5 } (34) 

ra = 41 (35) 
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For completeness, we note that the divergent terms satisfy 



CO 



T-(O) t(1) t(2) 

J n5 _|_ J n5 _j_ ^5 



(-1) 



55 



+ 



J, 



(0) 



55 



J, 



(1) 



55 



where we made use of 



L 3 



L 3 



- Of _ jk m —K^ - ik ( T (1) - T^ m ) 

2 1 — ^1^55 J tf " 2 m« — n v 55 "'m J 



and |H]). 

Using and (1351) we find 
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( 3T mn ) 
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(y) 
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— Z- 
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fa') 



a; 



where the tensor A mn is given by 
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and the sum over / in {ST}) runs from to 3. We have included the IZij terms together 
with the trace anomaly term (the last term on the right hand side of (j40p ) since the former 
will also induce a lne/L divergence in the stress tensor which is indicative of a conformal 
anomaly [M]. As observed before, we did not use a small k expansion to derive (139|) : in this 
sense it is non-perturbative. 

Suppose one is working in a regime where linearized hydrodynamics applies (i.e. small 
enough k). The first term in (1391) then relates to sound modes. The pole structure l/(k 2 — 
3a; 2 ) gives rise to a Mach cone in real space if the source is moving supersonically. The 
second term in (1391) . with the 1/u pole structure, is associated with diffusion of energy in 
the plasma, and in a suitable setup it signals the formation of a laminar wake far behind 
the source. It is commonplace for these poles to be shifted to slightly imaginary values by 
viscous effects, see equations f l63|) . Either the lo = ±y/3\k\ poles or the to = pole could be 
canceled by an appropriate numerator at leading order. If this happens, it corresponds to 
suppression of the sonic boom or the diffusion wake. 

In the near field, where the hydrodynamic approximation is not valid, we find that the 
Hu}\ completely cancel the pole structure in (1391) and therefore are not such good 

variables. In fact, in [351 [T4l [T6] it was shown how one may choose combinations of the 
components of the metric fluctuations whose equations of motion are completely decoupled 
from each other. While this is convenient when studying the solutions to the equations of 
motion at all scales, we find that the current choice shows the relation to hydrodynamics in 
a more transparent way. 

4 Large distance asymptotics of (ST mn ) 

As we have already remarked, if linearized hydrodynamics applies, {ST^) is expected to be 
subleading compared to other components of {ST mn ) by a single factor of k, because {ST^) 
has to do with viscous effects. If we set H^f\ = in fl39l) then at leading order in k the 
stress tensor is completely determined by the non-conservation equation (jH]) together with 
tracelessness (assuming no anomaly term A mn ). As a result, the constitutive relations of 
inviscid linearized hydrodynamics hold, likewise at leading order in k. One might think that 
the next correction in small k comes entirely from improving the constitutive relations by 
adding shear viscosity. What this section aims to show is that this is not quite right: in 
general, there is a non-hydrodynamical correction term J- mn , as described in (JH]). It enters 
at the same order as viscous effects if J(y) is of the same order in k as the four-force density 
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t(3) 
J m5- 

To find H^l, we need to consider some of the second order equations of motion, which 
appear in appendix lAl 



d vH {l]) + [— ) (1 - e ) i% } = -^J ( y) (1 - e y ) e 



.,2 



- ^k {i H j)0 (1 - e y ) 3/2 + 0(A; 2 tf mri ) (42a) 
o 

d 2 Y H 0l - dyH 0i = - Z l J w (l - e y )" 3 / 2 e y + 0{k 2 i H mn ) 

+ 0(ukiH mn ) , (42b) 

where we have defined a new variable 

F = In /i . (43) 

In the F variable, the boundary is located at Y = and the horizon at Y = —oo. The 
notation 0(kfH mn ) in (j42|) means any combination of two factors of spatial components ki 
with one factor of a component of the metric perturbation: for example, fcifei^oo would 
be such a term. Similarly, 0(ukiH mn ) in (142b|) means any combination of a factor of u, a 
component of the spatial momentum ki, and a component of the metric perturbation. In 
contrast, by O(k) we mean a term containing one factor of uj or fcj. 

We claim that in order to find (ST mn ) to an accuracy of one factor of k beyond the inviscid 
linearized hydro approximation, precisely the terms shown in (I42p need to be retained. This 
may seem counter-intuitive given that we have dropped factors quadratic in k^ from (|42|) 
while keeping factors quadratic in u. We will see in section section |4~31 that this seemingly 
uneven retention of 0(k 2 H mn ) terms is exactly what is required, because of the structure of 
the event horizon. 



4.1 Inviscid hydrodynamics and the non-hydrodynamical correc- 
tion 

Equations (I39p and (j42p . together with the assumption that the plasma is at rest at infinity, 
teach us that H mn must be 0(k~ n J fll/ ) with n > — 1 (we must allow n — — 1 due to the con- 
straint equations fl32|) ). Thus, to order O(k J IJ , u ), one may neglect 0{uj 2 Hij) and O(ukiH j) 
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terms, and (142 ap becomes 





X\-3/2 p Y 



(44) 



whose solution is 



H (ij) = Qij Y + 



fdy f dy'J {ij) {y'){l-e«)-^e 

JY J-oo 



(45) 



Now let's ask what happens when we plug the small z series expansions ([7]) and (|25|) into 
( 145]) and compare term by term. By comparing up to 0(z A log z) in (l4"5j) . one can recover the 
terms in (I26p that do not involve the Ricci scalar or tensor (more precisely, one can recover 
the mn = (ij) components of those terms.) These divergent terms are dealt with in section [3] 
and appendix O The 0(z A ) term in fj4"5l) gives us 



where Tij is defined in ( flTT) . 

To determine the integration constant Qij, we need to impose boundary conditions at 
the horizon which permit infalling waves but not outgoing ones. In the notation we're using, 



those which oscillate like e + ~i twY . This will be carried out and explained in more detail in 
section 14731 As a lowest order approximation (valid formally to order 0(k~ l J m b)) permitting 
only infalling modes implies that Hu^ asymptotes to a constant. Thus, according to (145]) . 
at this order = which, together with (I39p . gives us inviscid hydro up to the Tij terms 
and the trace anomaly term A mn . Corrections to at higher orders in small k will lead to 
viscous hydrodynamics. These will be dealt with in section H~3l 

4.2 Small k properties of J^(u,k,z) 



Before proceeding with an analysis of the near horizon boundary conditions, we make an 
aside to discuss in more detail what we mean by a small k expansion. Let's introduce the 
following notation for expansions in small k and small u: 




(46) 



this implies that we need to retain terms which oscillate like e * 



~ 1 - f iuiY but not 




oo 




(47) 
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where 

J$(\u,\k)=\*J$(u;,k). (48) 

In (1471) and (148]) we have suppressed the dependence on the radial coordinate z. In some 
instances will have delta-function support in Fourier space, meaning that uj will be 
localized around some u>o(k). This may happen when J^t) is not localized in time. To keep 
the discussion general, we consider a source whose stress energy tensor can be factored into 
an analytic and a non-analytic term: 

J^iuj, k, z) = A(u, k)j^(u, k, z) , (49) 

where A(u>, k) contains any non-analyticity. If J^iio, k, z) are analytic in u and k, we may 
simply set A = 1. The assumption that the probe source is localized in space amounts 
to having a smooth j^u, k, z) on the right hand side of (T451 : that is, j^ u (u,k,z) can be 
expanded in a Taylor series around k = 0. Actually, the definition of as presented in (|49|) 
is somewhat ambiguous since we may always rescale A and by appropriate factors of the 
momentum, for example, we may take A — > A/uj and — > u)j^ v . To fix this ambiguity 
we require that is of order k n with the smallest possible n. As noted in section [2] one 
may easily generalize these results to the case where J^ u is a sum of terms, each of the form 
shown on the right hand side of (T49|) . 

The superscripts n and n will be used fairly often in the following sections, not only for 
but for other functions as well. Thus, /( n > m ) specifies the n'th component of f(u,k,z) 
in a series expansion near the AdS boundary, and the m'th component of /^(to, k, z) in a 
small momentum/ frequency expansion. 

Shifting our attention to the near horizon asymptotics of J mn , we require that the stress- 
energy tensor of the source is causal in the sense that the (linear) response of the metric 
to the source will not be forced to have modes outgoing from the horizon, and that it is 
non- divergent. 

4.3 One order beyond inviscid hydrodynamics 

If the leading behavior of J^l(u, k) is of the same order or lower order (in a small momentum 
expansion) than the leading behavior of Fuj), then we can go one order beyond inviscid 
hydrodynamics and obtain the first viscous correction to it. For simplicity we assume that 
the leading behavior of is of order k°, as is the leading behavior of JF^/A. Our results 

(3) 

may be easily generalized to scenarios where starts at an order k n greater than that of 



1(3 



Tijj IS.. We will comment on this generalization where appropriate. Our strategy to obtain 
the inviscid corrections is to compare a small momentum/ frequency expansion of h mn (z), 
where H mn (z) = Ah mn (z), to a near horizon expansion of the solution. 

A small momentum/frequency expansion of huj\ can be read off from (145 p : 



z 2 f° f y 
lift) = qijY + -j j dy J dy'j (ij) (y')(l 



e'V'V, (50) 



where we have used Qij = qijA. To find g^, we consider the near horizon asymptotics of 
( )42|) (with H mn replaced by h mn and J mn replaced with j mn ). It may be intuitively helpful 
to think of (142j) as describing coupled oscillators evolving in a "time" parametrized by the 
coordinate Y (although in the AdSs-Schwarzschild geometry, Y is a spacelike coordinate) 
and perturbed by "forces" corresponding to the J mn terms. Near the horizon, the coefficient 
of huj\ in (|42aj) is constant, so huj\ "oscillates" with frequency Zqu/A. On other other hand, 
the hoi equations have no "restoring force" term, and the sign of the dyhoi term is such that 
the "motion" of h Qi is damped as one approaches the horizon. Thus the hoi tend to a finite 
limit as Y — > —00: 

lim h Qi (uj, k, Y) = -Wj(u, k) (51) 



for some Wj(u,k). The 0(k}h mn ) terms that were dropped from fl42b[) do not affect the 
conclusion (jBTj) because they are suppressed by a factor of e y , and because all the h mn must 



be assumed to be bounded as Y — > —00 in order for the linearized analysis to be valid. 
Using (|5"T|) the near horizon behavior of hnj\ takes the form 

h m = U^e-i**** + V tJ e^ Y + + ( ter S\hf L C oriIon iSh ) + °W ' ( 52 ) 

where Uij and Vij are integration constants and we have assumed that J(ij)e Y , and so also 
j(ij)e y , vanish at the horizon. Standard horizon boundary conditions are to allow only 
inf ailing modes: thus = 0. 

Comparing the small momentum limit of the near horizon asymptotics ( J52l) with the near 
horizon asymptotics of the small momentum limit (1501) implies that 

h m = U l3 (\ - \iz»uY^ + 2 -^p± + 0(k) = q l3 Y + (Fini^ o^caons) + (k) . (53) 



In ( 153|) we made use of the method of matched asymptotic expansions: we start with ( 1501) . 
which is valid for momenta small enough that the oscillating term at the horizon can be 
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Taylor expanded, — <C Y, and (1521) which is valid near the horizon, Y <C —1. The 
matching region is the intersection of these two intervals: 

<y<-i. (54) 



Provided zquj <C 1, which is to say u <C T, the region (154J) is non-empty. By matching the 
constant (F-independent) terms in ( |53l) we find 

jtF 1 ) = ^gg (55) 

By matching the terms proportional to y in (1531) we find 

^ = (56) 

q® = -\iz^U^ = ~iz k {i WP , (57) 



where in the second step we used fl55|) . 

To determine Wj , we consider the small momentum/frequency expansion of the equa- 
tion of motion for /ioi, ( I42bl) with if mn replaced with /i mn and J mn replaced with j mn . From 
( 1391) . the leading small momentum contribution to Sj, defined through the Poynting vector 
Si A = Si — — (T i), is of order 0{k~ l ). Thus, 

h Qi = sf l \e Y -1) + O(k ), (58) 

which gives us 

Wj = + O(k ) . (59) 



The traceless space-space components of the stress tensor in the absence of a trace anomaly, 
Hfy, can be read off of (1391, (155} and (150]): 

1 



= -^-^^(iSj) + Fij + (corrections) , (60) 

where we have used z = 1/nT. 

Plugging ( 1601) into ( 1391) we obtain ([9]), up to the aforementioned anomaly terms. The ex- 
plicit values for the energy density e = A ^e'-" 1 - 1 + e®^ and Poynting vector S = A ^s^" 1 ) + 
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are 



e (-D A = _ 3 ikdil ^05 0) A (61a) 

k 2 -3lu 2 

e (o) A = 3 ik ijS ,1)A ~ ^Jo5 1)a + fc^g ~ ik 2 k j s ( -~ 1) A/3nT 

k 2 -3u 2 

— ih i m -ik i m k luj 

FD A = + tifc ^ g '^j5 A (61c) 

W fc2 _ 3^2 

, m jj 5 3,T) A - iJL-J^ - fc i ife i sf i) A/127iT - fc 2 spW 47rT 
s- ; A = % 

UJ 

ikalpA - (ihjjf^A + kikj^ij - itfkjsf^A/ZitT) k { /u 

+ ^ =; L (61d) 

k 2 -3uj 2 

(3) 

Several remarks are in order. First, if the leading order behavior oij^ is of order k n then the 



superscripts (—1) and (0) should be replaced with (n — 1) and (n), where it is understood 
that s^" 1 ) = and e^ m ) = for m < n — 2. Also, if the leading behavior of Tij is of higher 
order in k than that of J m s then its contribution to ( 16TT) should be dropped. Second, we note 
that to leading order, a wake will exist as long as A does not cancel the 1/u pole of s|~ 1 . 
The expressions in ( 161bl) and (161d[) include terms behaving as 1/u 2 and l/{k 2 — 3uj 2 ) 2 



(coming, for example, from substituting the explicit expression for Sj into (I61bl) ). These 
terms can be understood as shifting the position of the sound and diffusion poles to slightly 
imaginary frequency, indicating viscous attenuation. Retaining O(k) accuracy and using 

Fm = J-m5 ~ 3~nm i (62) 



we find the following resummations of (ISTj) : 



UJ 2 — |/C 2 + 7^j;ik 2 0J 



O(kJ^) (63a) 



-^^F hydro + luhkiF^/k 2 K hydro - kikjF^Vk 2 
S = 3 ,° TT3 ■ 1 .... i i, + 0{kJ, v ) , (63b) 



UJ 



Ik 2 + ^fik 2 uj iuj - ^k 2 



which are precisely the expressions for the energy density and Poynting vector for a conformal 
fluid with rj/s = l/4n in the linear hydrodynamic approximation, sourced by an effective 
hydrodynamic four- force density F hydro . The effective four force density F hydro deviates 
from the actual four force J„J by the k n J r nrn term on the right hand side of fl62l) . It would 
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be interesting to see how higher order corrections arise [29j 12H] • Obtaining (|63|) completes 
our demonstration of in the case where J^l is of the same or lower order in a small k 
expansion than J 7 ^. 

From our construction, we can see how the constitutive relations arise from the Einstein 
equations. Near the horizon, the traceless space-space components of the metric are coupled 
to the Oi components of the metric in such a way that imposing no outgoing modes from the 
horizon induces the constitutive relations among the remaining integration constants. 



4.4 Configurations without a wake 



In contrast to the discussion in section 14731 let us consider the case where the leading behavior 
of J-'ij at small k dominates over the leading behavior of J^. Then (H6j) and the discussion 
following it implies that 



(3). 
m5 1 



(64) 



The 0(Jj n l) corrections now become more difficult to evaluate since the near horizon behavior 
of is hard to access in this setting. At this point, our perturbative analysis allows us 
to study only the inviscid limit of such configurations. If the leading behavior of J-'ij/ A is 
of order k° then the leading terms in an asymptotic expansion of the energy density and 

^ to zero. Otherwise, 



Poynting vector around k = can be read off of (l6Tj) after setting j m5 



(3,n) 
m5 



0. 



if J~ij is 0(k n A), we need to replace (0) and (1) with (n) and (n + 1) and set j 
As in the previous section, we will assume for simplicity that J-'ij is of order k°A. The 
generalization to a higher order dependence on the momentum is straightforward. 

Let's look for the conditions for the absence of a wake, or a 1/u pole in the leading terms 
for the Poynting vector, once is subleading. If A doesn't introduce such a pole, then this 
could happen if k x Tij has an uj dependence which will cancel the 1/u poles in (I6T|) : using 
(fTSj) we can rewrite ki J-'ij in terms of integrals over the Jq^'s and over Jj, 



ikjJ 7 ^ 



-iuF 



d( [d c {J l5 hC 3 ) 
C 3 h 



ihF 



c 3 j . 



J i5 h(~ 



' =0 



(65) 
(66) 



Note that the last two terms cancel each other at order 0(k~ 1 J^l). This follows from 
expanding (I18ap in small k and solving for z. To emphasize this point we substitute 
with Ajfo' 1 ^. Next, we note that generically, Joih" 1 diverges near the horizon. Since J 7 is 
finite, this divergence must be canceled by a similar divergence which can only come from 
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the next to last term. Thus, we write 



ikjJ-'ij/A = —iuF 



u( 3 



Z{) 



If 



then the integral 



C 3 h 



(o) 

i 



- -ikiF 



lim j\l ] h = 



30i C 
h 



-AdC 



.3 ( 



"3 h 



+ jHt T) + 0{k»). (67) 



(61 



(69) 



must be finite for z > and the 1/u pole will be canceled by the u dependence of the first 
term on the right hand side of (161dl) . We end up with 



3lu 



e = 1- 



k 2 - 3co 2 



J, 



(3) 



A- 2 



05 



k 2 - 3cu 2 



Si = i- 



k i 



(3) 



J i — L -» 05 ' ~* 

k 2 — 3u 2 k 2 — 3cu 2 



n C 

c 



3uki 



k 2 - 3uj 2 



c* JOj 



Joih 1 



k 2 - 3uj 2 



Z ° dC J ft" 1 



(70) 



(71) 



Since the tu -1 pole has been canceled, condition (168]) is sufficient for the absence of a wake 
as long as A does not introduce such a pole. Comoving string configurations which do not 
reach the horizon and whose monopole structure of j^l vanishes (for instance mesons or 
baryonic configurations) fall into this category. 



5 Short distance asymptotics 

In the previous section we extracted the large distance asymptotics of Hmn by appealing to a 
small momentum approximation, kz <C 1. Here, we consider the short distance asymptotics 
of the solution, associated with large momentum: kz$ ^> 1. The kz Q — > 00 limit corresponds 
to taking the temperature of the black hole to zero; as long as we keep terms up to order 
0(k~ A ), we are essentially working in the AdS geometry whose line element is given by ([2]) 
but with h = 1. More explicitly, by choosing 

X = ( #(11), #(12), #(13), #(22), #(23), --T #00 + #*, #*, #01, #02, #03 J , (72) 
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the second order linearized Einstein equations (1241) may be written as 

4 



a~ 3 h- TH d z {X^K*) + VijXj = --Si (73) 



where 



n= (1,1, 1,1, 1,3/2, 1/2, 0,0,0) (74) 



S — ( J J(12)i J(13), J (22), J(23), ~ (^~faJ° + "^"J ' /T^ 00 ' ^ 01 ' ^ 02 ' ^ 03 



(75) 



and Vij is an 0(k 2 ) 10 x 10 matrix given in appendix [Al In the kzo ^> 1 approximation 
/i = 1 + 0(zq 4 ), so up to order (zok)~ 4 all the kinetic terms in (1731 have the same form and 
Vij is a constant matrix. The left hand side of (173]) reduces to the equations of motion for 
linearized gravity in an AdS background and the source terms on the right hand side of ([731 
need to be expanded up to order (zo^) _4 E 

To solve ([73]) in the kz ^> 1 approximation, consider the five eigenvectors, x\ °f 
whose eigenvalues are —k 2 = uj 2 — k 2 , described in appendix [B] From (]73|) . the combinations 

= XjXj satisfy the massless scalar field equation in empty AdS space 

(a~ 3 d z a 3 d z - k 2 ) Y { = -A X )Sj . (76) 



The Green's function for the operator on the left hand side of (I7b|) has been extensively 
studied in the literature, see for example [36J. Using the Green's function, the solution to 
([76]) may be written as an integral over the source term [37] . We follow a somewhat different 
path to solve (176"]) . which allows us to write the solution to the equations of motion as a sum 
instead of an integral. Let's assume that 

J^u, k, z) = A c (u, k)j^ v (u, k, z) . (77) 

where j^ u {uj, k, z) has an expansion in integer powers of u and k for large values of ui and k, 
whereas A c (u, k) can be arbitrarily non-analytic. We expand 

oo 

S i = £ S t )za + 0{k- 4 A c ) . (78) 

a=-l 



4 Higher order corrections coming from the source term may also be treated. This is discussed in [14]. See 
also |15j for an explicit application. 
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According to the discussion at the beginning of this section, we are then allowed to solve 
(175)) with Si replaced by the right hand side of (1751) and with h — 1. In appendix iBl we 
construct the solution W n (z) to the equation 

(a~ 3 d z a 3 d z - k 2 ) W n {z) = z n (79) 



with appropriate boundary conditions for any n > — 1. Thus, the solution to (ITbl) takes the 
form 

oo 

yi = - 4 *5 E W a (z)S$ a) + 0(k^A c ) . (80) 

a = -l 

Recall that we are not interested in the full solution to (1751) . but only in the fourth order 
coefficient of Y l : 

oo 

Y 1 ^ = -A X ) E W^S^ + 0(k~ 4 A c ) . (81) 

o=-l 

Using (H33D, (USD, (H36D and ([139]) we find 



^i 4) (^) 



n-2 



-£ + ^7* - = (82) 

k (-l)f 2 n ~ 5 r (|n + 1) T (|n - 1) otherwise, 



where 7^ is the Euler-Mascheroni constant. 

With Y 1 ^ in hand we can solve xixf = Y i{ - 4 \ the energy conservation condition (13*2|) . 
and the trace equation (|3T|) to get the energy momentum tensor in terms of J^l and Y l ^ 4 \ 
Defining 

00 

= -4 E ^ (4) J - + ^(^" 4 A C ), (83) 

a=-l 

and 

-^m = 4i5 — ^ ^mn (84) 

oo 

= jS + 4X>£ ) 1 (3-a)J# (85) 

a=0 

(where in the last line we made use of (j 18 aft , which is valid to all orders in a in the 0(k~ 4 ) 
approximation we are using) and making use of (1301) . the fluctuations of the stress energy 
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tensor take the form 

(ST mn ) = T mn + T mn + A. mn (86) 

with 

T m n = ^ ip'^m^n "I" Vmnk ) — -g^ (k m F^ + k n F^j + g^^iP (k m k n — Tj mn k ) . (87) 
and 

A mn — 24 ^2 ^ron + 24p {^Vmnk 2 + ^k m k n ) Jm" 1 ~ -gJ^l, + ^ (^m^n ~ Vmnk 2 ) (ST}) (88) 

where we have made use of ( !33l . The trace of -4. mn , (ST}), is given in (l3~5l) . Note that 

^ m ^L = F*. (89) 

6 Strings and point particle sources 

In this section we consider two specific sources: stringy sources whose motion may be deter- 
mined through the Nambu-Goto action 

Sng= I drday/^gC NG C NG = -- — ■ - J^g (90) 

J lira 

with g a p = d a X^dpX u G llv , and point particles with mass m > whose action is given by 

^particle = / dfj Z/ particlo -^particle = G^ V X^X 2^ ^' (^"0 

where e is a Lagrange multiplier. In both instances the four-force induced by these objects 
in the boundary theory can be calculated for generic configurations. We find that a drag 
force acts on the string endpoint if it is located on the asymptotically AdS boundary while 
point particles, which are localized on the boundary only at some time t = TW,\, supply an 
impulse to the plasma at that time. 

6.1 Strings 

Strings whose endpoints are on the asymptotically AdS boundary are dual to infinitely 
massive quarks [201 12], mesons [381 121], or baryons if one introduces D5-brane baryon 
vertices (321 HQ]. 
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The energy momentum tensor (J4]) following from the Nambu-Goto Action (190]) is given 

by 

L 3 f 

J mn (t,x,z) = -- — - / da^a~H {3 \x- X(t,a))5(z- Z(t,a))d a X^d a X u (92) 

— /I LX J 

where we have used a gauge where t = r so that A M = \t X 1 X 2 X 3 Zj ■ We'll be 
interested in the boundary four-force density generated by the string. Therefore, we focus 
on the near boundary behavior of J m s where we can, at least locally, use the gauge Z = a. 
A near boundary expansion gives us 

J 05 = ^arV^f - X(t, z ))^= + ... (93) 
lira v~9 

L 3 , 1 (gaoX i '(t,z)-g w X%zj) 

^ = —a^h-WXS - X(t, z))± = L + .... (94) 

LTxa' v~9 

Defining Pj as the a component of the world-sheet current, 



8S NG 1 a 



2 



^g 00 X i '-g 01 X^ (95) 



5X U 2vra / ^/=^ 
(which is not necessarily a constant) we find that near the boundary, 

^string = _ Z 3 S (3)^ _ p. + (gg) 

^ring = z 3 6 (3)^ _ ^ ^p. + (g7) 

Using ([6]) and assuming a finite Pj on the boundary, we find that the four force is localized 
at the string endpoint(s) X % = X^ and is given by 

J { 3, = Y, 5i3) @-Xw)^ (-v- p p) ( 98 ) 

b z-> 

where V is the velocity of the string endpoint, V = \im z ^oX, and the sum is over the 
endpoints of the string which reach the boundary. A trailing string will clearly induce a 
wake since the total drag force is non vanishing. The total drag force acting on the center of 
mass of mesonic and baryonic configurations (which were recently introduced in this context 
in [H]) vanishes and therefore, from (191!) and the condition (|68|) a wake will not form (at 
least to leading order.) 
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In order to obtain the J- mn corrections, we need to consider an explicit string config- 
uration. The large distance and short distance asymptotics of (5T mn ) for various string 
configurations has been studied in [TH H5l HZl HH [221 CHI HH1 HH [33] . We have checked that 
our main result (Q fits with the expressions in the literature. The reader is referred to these 
references for the detailed structure of the stress-energy tensor. 

In sections [3] and Appendix Owe argued that the stress-energy tensor contains divergent 
terms which are localized to the extent that the probe-source is localized. For the case of 
string-configurations we will now show explicitly that these divergences can be associated 
with the infinite mass of the quark dual to the string endpoint. Since Pj and the velocity Vi 
are finite at the asymptotically AdS boundary, we find from (I95p that near the boundary, in 
the Z = a gauge, X-(t, z) is at least of order z 2 . This implies that 

Soo = -« 2 (1 - V 2 ) + O(z ) (99) 
<7io = O(z ) (100) 
gn = a 2 + O(z ). (101) 

Thus, to leading order in z, we find that a string endpoint which reaches the AdS boundary 
will induce an O(z) contribution to J mn of the form 

JW = M.eVT^U m U n 5®(2-X {b) ) (102) 

where 



is the four velocity of the quark and 



■1 V^/Vl^V 2 (103) 



M s = (104) 



is the mass of a static quark dual to a string ending on a D7-flavor brane [12] a distance 
z = e from the boundary [20]. In the e — > limit, or in the absence of D7-branes, the mass 
of the quark becomes infinite. With ( 11021) . we find that the divergent contribution of the 
stringy source to the stress energy tensor takes the form: 

{T e mn ) = ^M 9 VT=Wu m U n 5<- 3 \$-X (b) ) (105) 

b 

which is precisely what we would expect from a quark of (infinite) mass M s moving with 
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a 4- velocity U m . Naturally, it is localized at the position of the quark. Our result f 1 1 5 j) is 
an extension of the one obtained in [TO] for the special case of a trailing string moving at 
constant velocity. In the case of a stringy source there are no other divergent contributions 
to the stress-energy tensor. 



6.2 Point particles 

Here we discuss in some detail pointlike particles in the AdS black hole background. These 
can be thought of as approximations to quantum states of closed strings propagating in the 
bulk, which represent color-singlet quasi-particle excitations of the gauge theory — heuristically, 
glueballs which are in the process of thermalizing with the medium. Before computing the 
stress tensor, let us introduce 

BL 1 

Pm = -3T- = ~GmuX V (106) 

dX m e 

which is conserved since is independent of X m . Thus, we are free to define 

P m = (-E P 1 P 2 P 3 ) (107) 

with E and Pi constants. The stress tensor of the point particle following from varying the 
action f[9~Tl) with the conventions we introduced in (Tj|]) is 

I 3 1 

jmn = ^ =S (3), X _ X(t))S(z ~ Z(t))-X»X V (108) 

V-G e 

where we have used the gauge X M = X 1 X 2 X 3 Zj ■ Using (11061) to evaluate the 
near boundary asymptotics of J m5 we find 

J m5 = z 3 6^(x - X(t))5(z - Z(t))P m Z + . . . (109) 
= z 3 5®(x - X(z))S(t - T{z))P m + ..., (110) 



so that 

r \( b ))o*->{x - VL( 6 )J-r m 



42 = 5{t - T {b) )5^\x - X (b) )P m (111) 



where Tn* and X^ are the time and place where the particle reaches the boundary. 

To study the stress tensor of the plasma in response to the glueball we need to restrict 
ourselves to a more specific particle configuration. Consider a massless particle traveling 
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along x 1 



x 



with E > Pi > 0. Using the conservation equations (11061) and the equation 



of motion for e which follows from the variation of (1911) . one finds that the trajectory of such 
a particle satisfies 



X 1 



h(Z) X /l-^h(Z) 



(112) 
(113) 



where we have chosen a solution representing a point particle moving from the boundary to 
the horizon. Setting X 1 (0) = Z(0) = we get 



j£l=6(t)5^(x)(-E ?! o). 



(114) 



According to our discussion in section 14.11 this is enough to determine the leading long 
distance components of the stress tensor and the sub leading inviscid corrections to it. The 
non-hydrodynamical subleading corrections manifest themselves in terms of J r rnn , defined in 
(ITT]). From flTUSD, (TO and (JTTJ]) we find 



_ P 2 Z 3 ll(z) iuJ T(z)+ik 1 X 1 (z) 



(- 2 

3 







o o\ 

o h 



;ii5) 



where T(Z(t)) = 1. It follows that 
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•J %; 



8E 



f- 2 - 

3 
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lie) 



The energy momentum tensor of this configuration can be read off of © by using (fT3"]) , 
( |63|) . (I114p and (11161) . The resulting expressions are somewhat long as are their counterparts 
in real space. To get a flavor for the dynamics of the decaying glueball, consider a simplified 
trajectory where Pi = 0: a massless particle located at the boundary at t — and moving 
straight down into ArfS'5-Schwarzschild with x 1 = x 2 = x 3 = 0. In this case, the leading 
term for the energy density reads 



e(cu, k) 



3iuE 



k 2 — 3u> 2 — 3iT s k 2 uj 



;ii7) 
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where in the denominator we have added viscous corrections, T s = which enable us to 
carry out the Fourier transform in the to coordinate via contour integration. When Fourier 
transforming in the momentum coordinates we keep k small but allow kt to be large [2]. We 
find 



e(t,x) 



E 



t 



(r+c 3 t) Z 



(r-c a ty 
2F a t — g 2F a t 



e(t). 



118) 



4?rr " y/W^i 

where O(t) is the step function which vanishes for t < and r 2 = x l x l . Up to the step 
function, equation (j!18p coincides with the kernel for the sound mode contribution to the 
energy in the linear hydrodynamic approximation which was calculated in [2]. This is not 
surprising since the point particle provides an instantaneous impulse to the plasma. 

For the leading short distance asymptotics of the stress tensor, we need to compute T mn 
defined in f l83|) . Using rj = \J E 2 — P 2 , we find that the leading contribution to jF^ n reads 



r 



c 

mn 



-4- 



E 



Pi 




/ E 

-Pi 



\-v o 



11 
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-v 

E 'I 





p2 



^(_l)«/22 a - f T Qa+ 1^ 
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■ LuE~k 1 P 1 



a— 3 



3)! 



;ii9) 



The sum on the right hand side can be evaluated explicitly: 
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mn 



(5 + 2Qe , Ssinh" 1 ^) 
8g(l + e 



2\2 



8g(l + e 2 ) 5/2 
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-Pi 








— T] ' 
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p2 



(120) 



where £ 



" g -^ lPl and k 2 

r]k 



-UJ 



k 2 . By plugging f)120p into fl86l) we obtain the leading 



short distance asymptotics of the stress tensor. The expressions we find are somewhat 
complicated, so we will not reproduce them here. Instead, we pass to a special case: we 
consider once again the Pi — > limit and assume that uo has a slightly positive imaginary 
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part. Then 

e(uj, k) = i- 

\k\ \ w 



e(u,k) = z^tanh" 1 ( M ) . (121) 



Fourier transforming (I12ip we find 



e(t,x) = ^8(t-r) (122) 



where we have used Im{u;} > 0. We note that this expression can be obtained in a rather 
straightforward way by appealing to the symmetries of the problem. Using energy conserva- 
tion, the 0(3) symmetry in the directions, causality, and the fact that the point-like 
particle moving in the AdS bulk perturbs only a light-like component of the metric [13], we 
find that 



( ST uu) = 7Z3*(«) (123) 



E 
47rr s 
with u = t — r. 

From (11221) and (11181) we see that the energy density of a glueball which is "injected" 
into the plasma at t — will initially propagate outward with the speed of light, but at 
large distances, where hydrodynamics kick in, it will propagate with the speed of sound. We 
note that (11221) is in some disagreement with a result obtained in |37j . where it was claimed 
that the boundary stress energy tensor related to a massive particle moving in empty AdS 
vanishes. 
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A Linearized Einstein equations in the AdSs-Schwarz- 
schild background 

The second order Einstein equations (I24p for the linearized metric fluctuations in an AdSs 
black hole background (1221 are given by 



- 3 h~ n >d z (X>a 3 h n >) + V i0 X 3 = --Si 



(124) 



where X iy rii and Si are defined in (1721 . fl?4"|) and (1751) and is an C(/c 2 ) 10 x 10 matrix 



v = v + v 1 + v 2 



(125) 
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B Linearized Einstein equations in empty AdSs 

The Einstein equations in AdSs can be read off of (I124p by setting h = 1. In this limit, the 
matrix V T becomes a constant matrix. We will be interested in its five linearly independent 
eigenvectors 
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with eigenvalues — k 2 = uj 2 — k 2 . The other five eigenvalues of V T are zero with geometric 
multiplicity 4. 

By defining Y l = XjXj, we obtain the massless scalar field equation 

a~ 3 d z (a 3 d z Y l ) - k 2 Y { = -^Sj. (130) 
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Our method of solving f 1 1 3 j) is to first solve 



a~ 3 d z {a 3 d z W n ) - k 2 W n = z n , (131) 
and then use the solutions W n to construct the solution to (11301) through 

oo 



a=-l 



with Sj = J2 a Sj a) z a ^ To solve ([m]), we first define W n (z) = z 2 (ik)- n w n (ikz), so that ffTBTD 
takes the canonical form 

C 2 <(C) + C<-(4-C>n = C (133) 

where ( = ikz. The solution to (11331) with n = 3 is w 3 (() = ^H 2 (C), where H 2 (() is the 
Struve function whose properties have been studied and are tabulated in the literature. See 
for example |44j . Also note that the solution to (11331) with n = 4 is = ( 2 and the solution 
to (11331) with n = is wo = (~ 2 . By induction, we find that for n > 3 a solution to (11331) is 
given by 

w 2m {() = Am(C) n = 2m (134a) 

w 2m -i(C) = P 2m -i(C) + (-l) m (2m - 5)!!(2m - 1)H|#2(C) n = 2m - 1 (134b) 

where P n (C) are the polynomials 

puo = E(- 1 ) fc4fc ? "mu /V 2m " a " 2 ( 134c ) 

£^ (m — «)! (m — k — 2)1 

P m-Vr n fc ( 2m ~ 1 ) !! (2m -5)!! 2m _ 2k _ 3 

JW-i(C) - l^(-l) (2m _i_2]fc)!!(2m-5-21fe)!! C ' (134d) 



5 Whcn we will discuss the empty AdS geometry as a low temperature limit of the black hole geometry 
then we'll need to consider only those terms in Sj which are of leading order in a large momentum expansion. 
See section [5] for details. 
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Similarly, for — 1 < n < 2 we find 

w 2 = l + 4C 2 (134e) 
wi = ~C + \h 2 (Q (134f) 
w = C 2 (134g) 

«'-i = ~r 1 -ic+ ^(0- (i34h) 

Equations (I134p are special solutions to the non homogeneous equation (I133p . The most 
general solution to (I133P are those given in (11341) in addition to the solutions to the ho- 
mogeneous equations J 2 (C) an d ^(C) which are Bessel functions of order 2. So we should 
take 

w n -> w n + a„J 2 (C) + b n Y 2 (C). (135) 

Since Y 2 = — ^C~ 2 + • • then the boundary condition H mn (0) = 0, which translates to 
lim^o C 2 ^n(C) = 0, implies that 

b = (136) 

b 2 = 7T (137) 

and b n = for all other n. The coefficients a n can be found by considering spacelike momenta, 
k 2 > 0. Then, 

H 2 (ikz) + aJ 2 (ikz) = —iL 2 (kz) — al 2 (kz) (138) 

where L 2 is a modified Struve function and I 2 is a modified Bessel function. Both L 2 and I 2 
diverge exponentially as we take z — >• oo, but the combination L 2 — I 2 remains finite. Thus, 

-i n = 2m — 1 
a n = { (139) 

n — 2m. 

Now that the solution to (11331) is available, we can easily solve fl 1 30[) . 

oo 

Yl = - 4 *5 J2 z 2 w a {zq)Sf\ikY a . (140) 

a=-l 

For completeness, we note that the remaining five equations of motion will leave us with 
five undetermined integration constants which, as described in section 13.31 are determined 
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by the five first order constraint equations. 



C Holographic renormalization 

The prescription for obtaining the boundary theory stress tensor (T mn ) from an asymptot- 
ically AdS bulk metric has been developed in [45j [M] following the basic prescription of 
[TT| [T2] and is given by 

(T mn ) = lim ; 2 - 9Sr ;: , (141) 

where g mn (x, z) is defined through the line element 

r2 

ds 2 = — (dz 2 + g mn (x, z)dx m dx n ) (142) 
z 2 

and z — > corresponds to the asymptotic boundary We take e to be small so that 

L 2 

lmn(x,e) = —g mn (x,e) (143) 
e z 

gives us the metric on a spacelike hypersurface of constant z = e close to the conformal 
boundary (and it is understood that 7 5m = 0.) The action S ren is given by 

Sren = SeH + Sj\j + SqH + S'ct- (144) 

The Einstein-Hilbert action Seh an d the matter action Sm are given in ([T]). The Gibbons- 
Hawking boundary term reads 



Sen 

^5 J z =e 



\ [ d^x^K (145) 

K 5 Jz=e 



with the extrinsic curvature, 



c 



and 



X M „ = = -^Iw (146) 



n» = ~5^. (147) 
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The boundary action, S ct , is constructed from the boundary values of the dynamical fields. In 
the probe approximation which we are working in, the only dynamical field is the boundary 
metric 7 mn so that S ct = S'ctM- More explicitly, 

Sct = ~^J z d S/7 (|a a + irLA fl ) (148) 

with Aa and Ar undetermined numbers and W is the Ricci Scalar associated with the 
boundary metric 7 mn . It is easy to convince oneself that, in the setup we're considering, any 
higher order terms in derivatives of 7 will not contribute to (T mn ). If a conformal anomaly 
is present then S c t can also depend explicitly on lne [34]. Usually a logarithmic correction 
to 7, nn at order 0(e 2 ) is indicative of such an anomaly. In our case, if j£l ^ or J® ^ 
such a logarithmic term will appear in the series expansion for H mn and we will find a trace 
anomaly for the boundary stress tensor (see (I39p .) However, since in our case the logarithmic 
terms are not induced by a deformation of the boundary theory, as in [M], we will need to 
deal with the In e divergences (and other divergences) in a somewhat different manner. 

C.l Evaluating (T%%) 

The program of holographic renormalization [4"6"l W7\ provides a prescription for choosing 
S ct such that all the correlation functions of the boundary theory will be finite. In the 
current context, this amounts to finding the two coefficients Aa and Ar in (11481) . Thus, it is 
much simpler to find the appropriate counterterm action S c t by requiring that the one-point 
function (11411) is finite in the e — > limit. This method has been used in [48] to find the 
counterterm action in the presence of irrelevant operators. Here, the benefit is that we do 
not have to deal with the unknown contribution of Sm to the action when evaluated on-shell: 

= £s jm wk^) {Seh + s - + Sch) (149) 

2 d 

= lim — - Sqh (150) 

-o^)9f n (l,e) 

L 2 

= - lim -5-^ (K mn - K'jmn) . (151) 

Using (11461) we find 



L 3 ( 3 g mn g ts d t g ts - d e g r 
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Rewriting the metric (1221) in the Graham- Fefferman coordinate system (11421) . 



9mn = Vmn + ^£^3 (-^mn)bath e + ^J^H mn + 0{e ) , (153) 

(11521) reduces to 

(T GH ) = lim ( —^H^ + J_d ( H _ „ + ( finite^) 

\ J mn/ ™ol K 2 e 4 ^^^V"™* ^^fi J + ^ terms/ • ^ i0 ^ 

In what follows we expand H mn in a power series, as in (1231) . but with an additional 
H mn (0) = Hmn 7^ term which we shall set to zero at the end of the day. The reason 
we temporarily allow such a term will become clear shortly. To express (11541) in terms of 
the energy^ momentum tensor of the probe source, we use (126]) to substitute the Hml's with 
the J, 



a ) 6 
in &U 



Working withe boundary metric 7 mn = ^2g mn instead of g mn , we find that the 



divergent terms in (T^f ) are given by 



12 ^fo'-^^^-L^'-iff^^ (155) 



where R? = j^R and i?^ n = R mn are the Ricci scalar and Ricci tensor corresponding to the 
boundary metric •y mn . 

C.2 Evaluating the counter term action 

Some of the divergences in (I155P can be removed by an appropriate choice of Aa and An in 
( 11481) . Comparing 

*3 ; h ^JL.A * = * ( ——j- + lAr (V - U^)) (ise) 



to (I155p . we find that to get rid of the leading divergences we must set 



A A = -3, A R = -i. (157) 



6 The astute reader might worry that equation (|26|1 has been constructed for the Hmn = case. As it is 
written, it is also valid for non vanishing Hmn (and thus non vanishing Rmli) as long as we arc working to 
linear order in H mn . 
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Had we've taken Hmn = 0, we would have concluded that 



A A = -3, A« = -l; 



(158) 




boundary metric, we should use the same S ct no matter what the boundary conditions on 
H mn (0) are (as long as they're finite.). Also, by setting J mn = we should get the coun- 
terterms needed for a pure gravity theory with a non flat boundary metric. These have 
been worked out in [31] and coincide with (11571) once we take into account the appropriate 
conventions. Thus, we must use (j!57p . 

At this point we are still left with a stress-energy tensor which includes power law and 
logarithmic divergences. In the pure gravity theory logarithmic divergences may be canceled 
by a boundary action containing four derivatives of the boundary metric 7. The reason these 
kinds of counterterms are available in the pure gravity case is that there, (R^n ~ \R 1 lmn) 
may be written as a certain combination of RZin , R 1 ^ and 7™ and this combination may 
be obtained by varying a higher derivative boundary action [31]. While we may add such 
higher derivative terms to (I148p . they will not cancel any of the remaining divergences in 
(11551) since the latter are independent of RZm\ R 1 ^ and 7mL Put differently, the probe- 
source acts as an external current in the equations of motion. As a result, the divergences 
it introduces can not be subtracted away by counterterms constructed only of the metric 
tensor. Luckily, since the Jmn's have finite spatial extent, these divergences can be thought 
of as divergent contact terms, and we expect that they are associated with formally divergent 
parameters of the source. 

After setting H mn (0) = and separating the divergent terms (T^ n ) from the finite terms 
of the stress tensor in (11411) . we find 



{T mn ) (2~mn) bath + (T mn ) + 7] mn H^ ^ + 



3L 3 



( fl(2) _ i R (2) 




(159) 
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with 



" f G ~ 7 - - ^ (2) ^«) ln ^) • ( 16 °) 
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